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Abstract: 

In  this  study  we  evaluated  the  efficiencies  of  various  thermomagnetic  harvesting  methods  and 
materials.  Initially  we  focused  on  producing  gadolinium  (Gd)  single  domain  nanostructures  which  we 
predicted  to  have  a  thermomagnetic  harvesting  efficiency  of  approximately  30%  of  the  Carnot  limit.  We 
successfully  deposited  epitaxial  thin  films  of  Gd  which  exhibited  a  spin-reorientation  response  and  saw  a 
reduction  in  Curie  temperature  when  compared  to  bulk  values.  Unfortunately,  we  were  unable  to 
produce  single  domain  elements  using  the  lift  off  methods  attempted.  This  was  attributed  to  poor 
sidewall  coating  of  the  capping  layer  which  led  to  a  non-uniform  Gd  morphology  and  stacking  faults  that 
significantly  reduced  the  saturation  magnetic  moment  of  the  material.  Upon  examining  the  other  single 
element  magnetic  materials  in  the  literature,  we  discovered  that  more  ordered  transitions  led  to  higher 
transduction  efficiencies.  Utilizing  this  phenomenon,  we  explored  the  effect  of  order-to-order  spin 
reorientation  harvesting  in  neodymium  cobalt  (NdCo5).  Analytical  modeling  showed  that  NdCo5 
operating  about  its  spin-reorientation  temperature  could  achieve  an  energy  density  greater  than  2 
MJ/m3  with  an  efficiency  of  approximately  22%  of  the  Carnot  limit.  We  attempted  to  produce  epitaxial 
NdCo5  by  the  same  sputtering  methodology  used  to  fabricate  the  Gd  films,  but  every  combination  of 
sputter  and  annealing  parameters  yielded  highly  textured  yet  non-aligned  crystal  morphologies.  While 
we  were  unsuccessful  in  producing  a  working  prototype,  we  garnered  significant  knowledge  about  the 
physics  governing  thermomagnetic  transformations  to  pave  the  way  for  future  work  in  this  area.  This 
report  overviews  the  body  of  literature  related  to  the  thermomagnetic  field,  provides  a  detailed  analysis 
of  the  UCLA  modeling  effort,  and  describes  fabrication  and  testing  techniques  used  to  create  and 
evaluate  nano-structured  thermomagnetic  materials. 


1.0  Background 


All  ferromagnetic  materials  experience  a  thermally  induced  second  order  phase  transformation  and 
transform  from  the  ferromagnetic  ordered  state,  where  the  magnetic  moments  prefer  to  align  with  one 
another  yielding  a  spontaneous  magnetization,  to  a  disordered  paramagnetic  state,  wherein  the  thermal 
energy  experienced  by  each  magnetic  spin  overcomes  the  aligning  force  and  the  material  no  longer 
appears  magnetic  to  an  outside  observer.  The  temperature  at  which  this  transition  occurs  is  known  as 
the  Curie  Point  (Tc).  By  oscillating  a  magnetic  material  about  this  temperature,  thermal  energy  can  be 
transduced  into  magnetic  energy,  where  it  can  be  harvested  as  useable  electrical  energy  utilizing  the 
principles  of  electromagnetic  generators. 

The  first  attempt  at  thermomagnetic  transduction  was  proposed  by  Nicola  Tesla1  in  his  1889  patent, 
Thermo-magnetic  Motor,  wherein  he  showed  how  a  thermally  oscillated  magnetic  system  could 
generate  mechanical  work.  However,  utilizing  such  a  system  to  generate  electrical  power  wasn't 
explored  until  Brillouin  and  Iskenderian2  in  1948  .  They  proposed  that  a  ferromagnetic  material  can 
perform  magnetic  work  on  a  coil  by  flux  oscillations  which  occur  as  the  material  is  thermally  cycled 
about  its  Curie  Point  (Tc).  Their  model  consisted  of  a  prototypical  Carnot  cycle:  heating  and  cooling  the 
device  in  the  presence  of  a  varying  magnetic  field,  wherein  electrical  energy  is  produced  by  the  flux 
change  in  the  magnetic  material  during  the  two  adiabatic  portions  of  the  cycle.  Their  phenomenological 
analysis  suggested  that,  in  the  limiting  case,  device  efficiencies  as  high  as  55%  of  the  Carnot  limit  are 
possible. 

One  of  the  major  obstacles  in  experimentally  verifying  this  theory  in  the  1940's  was  the  lack  of 
ferromagnetic  materials  whose  Tc  was  near  room  temperature.  Therefore,  this  area  remained  dormant 
for  the  next  decade.  In  1959,  Elliot  et  al3  experimentally  evaluated  gadolinium  (Gd)  for  use  in  a 
thermomagnetic  energy  harvesting  system  operating  at  room  temperature.  However,  when  actual 
power  density  calculations  were  made,  Elliot's  thermomagnetic  generator  failed  to  produce  any 
significant  power  outputs.  One  of  the  main  issues  cited  at  this  time  was  the  lack  of  sufficiently  strong 
magnets  to  bias  the  structures. 

Following  the  development  of  rare  earth  magnets,  in  the  1980's  Kirol  et  al4  proposed  a  regenerative 
thermomagnetic  cycle  for  increased  power  output.  This  regenerative  thermomagnetic  approach 
suggests  a  maximum  relative  efficiency  (r)Rei)  of  75  %.  However,  this  approach  was  never  experimentally 
attempted.  In  1991,  Solomon  et  al 5-7  designed  and  evaluated  a  thermomagnetic  generator  on  the  scale 


of  a  large  power  plant  (e.g.  Hoover  Dam)  using  Yittrium  Iron  (Y2Fei7)  as  the  ferromagnetic  switch.  Y2Fei7 
was  chosen  since  its  Curie  temperature  could  be  varied  by  compositional  modifications.  While  advancing 
the  field  of  regenerative  thermomagnetic  energy  harvesting,  the  negligible  power  output  limited  the 
usefulness  of  this  structure. 
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Figure  1:  Illustration  of  concept  using  layered  thin  film  along  with  thermal  and  magnetic  fields 

Following  this  body  of  literature,  UCLA  began  conducting  studies  on  macro  and  nanoscale 
thermomagnetic  energy  harvesting  structures,  postulating  that  significant  improvements  in  efficiency 
and  energy  density  could  be  made  at  reduced  length  scales.  The  method  being  used  to  harvest  these 
magnetic  changes  and  convert  the  energy  into  a  usable  electrical  form  is  by  harnessing  a  multiferroic 
composite  of  magnetic  and  piezoelectric  material.  This  system,  illustrated  in  Figure  1,  consists  of  a 
layered  beam  element  containing  a  ferromagnetic  switch  (e.g.  Gd,  which  has  a  room-temperature  phase 
transformation),  a  classic  ferroelectric  like  lead  zirconate  titanate  (PZT),  a  heat  source  containing  a  hard 
ferromagnetic  material  (e.g.  neodymium  iron  boron  (NdFeB)),  and  a  heat  sink.  At  the  heat  source,  the 
NdFeB  magnet  produces  a  magnetic  field  that  induces  an  attractive  force  in  the  Gd  material,  due  to 
magnetostatics  and  domain  wall  motion.  As  the  beam  bends,  the  beam  produces  charge  on  the 
piezoelectric  surface,  which  is  harvested  as  electrical  energy.  When  the  Gd  touches  the  hot  source,  it 
heats  up  and  undergoes  a  second  order  transformation  into  a  paramagnetic  (i.e.  non-magnetic)  state, 
releasing  the  beam  from  the  hot  source.  This  returns  the  beam  to  the  environmentally  available  cold 
sink  and  generates  additional  charge  on  the  piezoelectric.  Once  in  contact  with  the  cold  sink,  it 
undergoes  a  reverse  second  order  transformation  back  into  a  ferromagnetic  state  that  is  once  again 
magnetically  attracted  to  the  environmentally  available  hot  source.  This  cyclical  process  repeats  itself 
and  produces  electrical  energy  output.  UCLA  has  demonstrated  the  concept  at  the  mesoscopic  level  (1 
cm);  however,  the  measured  efficiencies  of  the  current  device  are  on  the  order  of  10-20%  of  Carnot,  i.e. 
similar  to  the  Seebeck  thermoelectric  device. 
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Figure  2:  Illustration  of  domain  walls  and  calculation  of  efficiency. 

The  UCLA  device  operation  is  represented  by  exchanging  thermal  energy  for  magnetic  energy  as 
illustrated  in  Figure  2.  Heating  the  ferromagnetic  material  by  AT  (i.e.  through  its  Curie  temperature) 
requires  a  specific  amount  of  thermal  energy  as  illustrated  in  Figure  2.  This  represents  the  input  energy 
into  the  device.  When  the  ferromagnetic  material  is  cooled  below  its  Curie  temperature,  the 
ferromagnetic  material  stores  some  of  this  input  thermal  energy  as  magnetization  in  the  material.  The 
remainder  of  the  thermal  energy  is  delivered  to  the  cold  source  as  waste  heat  and  represents  thermal 
energy  loss.  The  magnetic  energy  stored  in  the  ferromagnetic  material  represents  the  available  energy 
to  harvest  as  electricity. 


Figure  3:  Isothermal  magnetization  curves  and  thermomagnetic  cycle. 


An  idealized  thermomagnetic  cycle  which  utilizes  UCLA's  multiferroic  device  design  is  shown  in  Figure  3 
along  with  a  plot  of  magnetization  as  a  function  of  magnetic  field  (M-H  plot)  for  two  different 
temperatures  straddling  the  Curie  Point.  The  four  illustrations  shown  in  Figure  1  contain  a  magnetic  heat 
source  which  generates  a  magnetic  field  to  attract  the  magnetic  working  element  and  a  cold  heat  sink  to 
shuttle  packets  of  thermal  energy  to.  The  process  can  be  described  as  follows:  Initially,  at  point  1,  the 
temperature  of  the  material  is  above  Tc  (at  Thot)  and  is  therefore  non-ferromagnetic.  By  bringing  the 
material  into  thermal  contact  with  a  cold  reservoir,  the  temperature  drops  to  TC0|d  (TC0|d  <  Tc)  and  the 
material  becomes  ferromagnetic  at  point  2.  During  this  process,  the  material  is  subjected  to  a  magnetic 
field  which  increases  the  magnetization  of  the  material  to  point  3  on  the  M-H  plot.  Simultaneously,  the 
material  is  brought  into  thermal  contact  with  the  hot  reservoir  where  the  temperature  again  increases 
to  Thot  at  point  4  and  the  material  becomes  non-ferromagnetic,  returning  to  the  cold  side  at  point  1  and 
thus  completing  the  cycle. 

In  order  to  compute  the  amount  of  energy  available  per  cycle  we  needed  to  consider  the 
thermodynamics  of  the  system.  The  thermomagnetic  cycle  involved  has  a  net  work  which  can  be 
defined  as  the  area  between  isothermal  magnetization  curves.  If  we  consider  the  following  work 
expression  (expressed  in  cgs,  where  the  units  are  erg/cm3) 

Wnut - —  f  HdB - -if  (HdH  +  AjtHdM  IT,!!))  (1.1) 

Ajt *  cycle  4 JT*  cycle '  \  /  / 

Where  B  is  the  magnetic  induction,  H  is  the  magnetic  field,  and  M  is  the  magnetization.  Since  the 
integral  of  the  HdH  term  is  zero  over  a  complete  cycle,  Equation  1.1  reduces  to 

w--y^HdMdM).  (i.2) 

By  utilizing  energy  conservation,  the  heat  that  must  be  input  into  the  system  is  defined  as  follows: 

THot  THot 

ft,  -  a,  +  -PfCp  (T)dT  +  TfdSm  +  Wout  .pfCr(T)dT  (1.3) 

Tcold  Tcold 

Where  p  is  the  density,  Cp  is  the  specific  heat  at  constant  pressure,  and  Sm  is  the  magnetic  entropy 
associated  with  irreversible  domain  wall  motion  and  nucleation.  For  small  H  fields,  the  integral  of  the 
specific  heat  term  is  much  larger  than  the  other  two  terms  relative  to  the  magnetization  of  the  material, 
and  thus,  the  heat  equation  can  be  simplified  as  prescribed  in  equation  1.3.  With  the  work  out  and  heat 
in  equations  defined  for  the  system,  the  absolute  efficiency  can  be  described  as  follows: 


(1.4) 
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The  absolute  efficiency  can  be  compared  to  the  limiting  case  of  the  Carnot  efficiency  (ncarnot)  operating 
between  the  bounding  temperatures,  where  the  Carnot  efficiency  represents  the  maximal  efficiency  of  a 
heat  engine,  by  defining  a  relative  efficiency. 
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Therefore,  the  efficiency  of  the  system  is  the  amount  of  magnetic  energy  stored  in  the  ferromagnetic 
material  divided  by  the  amount  of  thermal  energy  required  to  raise  it  by  AT  (i.e.  AM/AT).  Maximizing 
the  device  efficiency  requires  maximizing  the  magnetic  energy  in  the  ferromagnetic  material  while 
minimizing  the  amount  of  thermal  energy  input  required  to  induce  a  phase  transformation  (i.e.  minimize 
AT). 

2.0  Harvesting  About  the  Curie  Point:  Single  Domains 

After  a  considerable  amount  of  research,  UCLA  believes  the  critical  feature  limiting  the  efficiency  is  the 
domain  wall.  The  premise  of  this  claim  is  that  as  the  magnetic  material  is  reduced  beyond  a  critical 
volume  that  it  becomes  energetically  favorable  for  every  magnetic  spin  to  be  perfectly  aligned,  even 
after  the  initial  biasing  field  is  removed.  This  state  is  known  as  Single  Domain  (SD)  and  can  be  further 
characterized  by  a  lack  of  domain  walls  and  a  remenant  magnetization  of  greater  than  90%  of  the 
materials  saturation  magnetization.  We  initially  focused  on  Gd,  due  to  its  room  temperature  Tc,  and 
evaluated  the  magnetic  energy  output  change  as  a  function  of  thermal  input  for  different  sizes  with  a 
focus  on  single  domain  limit. 
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Figure  4:  Isothermal  magnetization  curves  of  gadolinium  at  an  applied  filed  Ha=3000  Oe 

Figure  4,  Hsu  et  al8,  shows  Isothermal  magnetization  curves  for  Gd  in  its  polycrystalline  film  and  SD 
forms.  It  can  be  seen  that  the  area  between  curves  represents  the  work  produced  during  a  single 
thermomagnetic  cycle  operating;  in  this  case,  the  temperature  difference  between  the  curves  was 
chosen  to  be  AT=5  K.  At  a  length  scale  below  the  SD  limit,  the  remenant  magnetization  approaches  the 
saturation  magnetization  at  3000  Oe,  and  the  work  potential  of  the  system  increases  dramatically  from 
24.2  kJ/m3  to  69.2  kJ/m3. 

We  began  our  experimental  work  by  attempting  to  fabricate  single  domain  elements  of  Gd,  as  this 
element  could  be  characterized  about  both  its  Tc  and  SR  regions.  This  was  initiated  by  attempting  to 
evaporate  Gd  using  the  CHA  Industries  Solution  Evaporator  available  at  UCLA's  Integrated  Systems 
Nanofabrication  Cleanroom  (ISNM).  Even  though  these  samples  were  capped  with  about  10  nm  of  self- 
passivating  platinum  (Pt),  we  were  unable  achieve  the  magnetization  values  present  in  the  literature. 
We  determined  that  this  was  due  to  the  fact  that  the  cryopump  used  to  maintain  the  vacuum  was  only 
rated  for  5*10’7  Torr,  which  was  insufficient  to  prevent  oxidation  from  occurring  in  Gd  during  the 
evaporation  process.  As  evaporation  systems  are  also  unable  to  produce  crystalline  films,  we  changed  to 
an  ultra  high  vacuum  (UHV)  sputtering  system.  The  Denton  UHV  sputtering  system  was  capable  of 
achieving  5*10"8  Torr  and  was  run  in  DC  sputtering  mode  in  an  Argon  (Ar)  atmosphere  with  a  pressure  of 


10‘3  Torr  with  a  flow  rate  of  20  SCCM.  The  99.9%  pure  Gd  target  was  sputtered  at  30  W  and  a  deposition 
rate  of  ~1  A/s  to  a  thickness  of  25  nm  and  was  then  capped  with  20  nm  of  Pt  before  breaking  vacuum. 


Figure  5:  Left  -  XRD  results  for  20  nm  sputtered  Gd  thin-film  (a)  and  the  Silicon  substrate  (b).  Right  -  ZFC-FC  results  for  a 
temperature  range  of  10  to  350  K. 

The  sputtered  film  was  examined  by  x-ray  diffraction  (XRD)  to  test  for  crystallinity  and  super  conducting 
quantum  interference  device  (SQUID)  magnetometry  to  examine  how  the  magnetization  changed  with 
temperature  using  the  zero-field  cooling/field  cooling  (ZFC-FC)  technique.  These  results  are  presented 
in  Figure  5.  It  can  be  seen  from  the  XRD  intensity  plot  that  the  Gd  film  showed  a  high  degree  of 
crystallinity  with  a  sharp  peak  in  the  [002]  direction,  corresponding  to  the  C-axis  direction  of  the 
hexagonal  crystal.  The  ZFC-FC  measurement  was  taken  in  the  in-plane  orientation  using  a  SQUID  and 
showed  a  distinct  hump  between  220  and  260  K,  indicating  that  the  magnetic  moment  rotated  into  the 
direction  of  measurement.  This  anomaly  is  consistent  with  measurements  for  single  crystal  Gd9'10  and 
shows  that  the  film  was  deposited  with  sufficient  epitaxy  to  achieve  spin  reorientation,  which  is  an 
interesting  distinction  between  this  work  and  the  previously  published  work  by  Berger  et  al.10  wherein 
no  temperature  dependent  magnetic  rotation  occurred.  Their  explanation  was  that  sufficiently  thin  films 
exhibit  such  a  large  shape  anisotropy  that  rotation  into  the  out  of  plane  direction  should  be  entirely 
energetically  unfavorable.  We  believe  that  grain  boundaries  observed  in  their  film  acted  as  effective 
defects  which  interfered  with  the  long  range  order  required  to  produce  spin  reorientation,  but  that  the 
large  degree  of  epitaxy  observed  in  our  films  facilitated  the  spin  reorientation  response. 
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Figure  6:  Arrott  plot  for  25  nm  sputtered  Gd  thin  film 

Figure  6  presents  the  Arrott  plot  consisting  of  M2  vs  H/H  values  obtained  from  the  isothermal 
magnetization  data  of  our  epitaxial  Gd  thin  films.  The  curves  are  generated  from  measurements 
conducted  at  five  different  temperatures.  The  Curie  temperature  is  canonically  defined  from  an  Arrott 
plot  as  the  temperate  whose  curve,  when  extrapolated,  intersects  the  origin.  For  this  Gd  film,  it  was 
found  that  the  saturation  magnetizations  measured  at  the  lower  temperatures  are  similar  to  those 
presented  by  Dan'Kov  et  al.12;  however,  the  measured  values  became  disparate  at  higher  temperatures. 
The  Arrott  plot  from  our  film  shows  that  the  Curie  temperature  of  this  film  is  approximately  278  K, 
representing  a  roughly  15  K  shift  from  the  bulk  critical  point  at  ~298  K,  which  we  verified  from  our 
evaporation  source  ingot.  In  the  literature,  shifts  in  the  Curie  temperature  were  also  observed  in  various 
epitaxial12  and  sputtered13  Gd  films.  We  attribute  the  change  of  Tc  in  our  sample  to  internal  stresses 
which  were  introduced  during  the  sputtering  process.  This  process  of  shifting  Tc  is  an  interesting  one 
because  it  shows  that  by  modifying  the  crystal  and  magnetic  ordering  of  a  material,  in  this  case  by 
depositing  a  very  thin  film  epitaxially,  it  is  possible  to  tune  the  thermomagnetic  behavior. 


Figure  7:  Left  -  Isothermal  magnetization  curves  of  Gd  nanobar  array.  Right  -  (a)  TEM  image  of  Gd  nanobar,  (b)  Diffraction 
patter  of  nanobar  cross-section,  (c)  HRTEM  of  nanobar,  (d)  Select  area  FFT  image  with  FCC  zone,  (e)  HRTEM  image  of 
nanobar,  (f)  Selected  area  FFT  image  of  HCP  zone. 

When  the  same  deposition  process  was  used  in  combination  with  a  lift-off  procedure  to  produce  nano¬ 
structures,  in  an  attempt  to  create  single  domain  elements,  the  saturation  magnetization  and  coercivity 
results  proved  unsatisfactory.  These  results  can  be  seen  on  the  left  side  of  Figure  7.  From  the  M-H 
curves  it  was  observed  that  the  saturation  at  53  K  is  less  than  500  emu/cm3,  when  it  should  be  closer  to 
2000  emu/cm3,  and  has  a  coercivity  of  less  than  500  Oe.  This  most  probably  indicates  a  non-uniform 
crystallinity.  It  was  also  expected  that,  due  to  the  magnetocrystalline  anisotropy  change  across  the  spin 
reorientation  transition,  there  would  be  a  change  from  easy  to  hard  axis  as  the  temperature  was 
increased,  however  that  was  not  observed  in  the  nanostructres.  Given  these  disparities,  we  examined 
the  nanobars  under  a  high  resolution  transmission  electron  microscope  (HRTEM),  Figure  7  Right,  where 
we  noticed  that  the  Pt  capping  layer  failed  to  coat  the  sidewalls  after  lift-off,  leading  to  a  dramatic 
reduction  in  the  amount  of  ferromagnetic  material  present  by  oxidation,  and  that  different  regions  of 
the  bars  had  different  crystal  structures;  the  cubic  (FCC)  zone  seen  in  subset  (d)  produces  significantly 
less  magnetocrystalline  anisotropy  than  the  hexagonal  (HCP)  zone  in  subset  (f). 


The  reduction  in  saturation  magnetization  along  with  the  absence  of  shape  anisotropy  induced  single 
domains  in  the  Gd  nanobars  can  be  attributed  to  oxidation  as  well  as  non-uniform  morphologies  in  the 
films  crystallinity.  The  presented  results  show  that  the  crystallographic  ordering,  and  uniformity  of  the 
film  has  a  substantial  impact  on  the  magnetic  anisotropy  and  magnetic  moment  of  Gd  nanostructures. 

We  also  investigated  a  wide  range  of  single  element  materials  where  experimental  data  was  readily 
available.  Table  1  outlines  some  of  the  materials  explored  in  their  single  and  polycrystalline  states,  their 
Curie  temperatures  (Tc),  as  well  as  the  calculated  relative  efficiencies  when  operating  about  Tc. 


Table  1:  Thermomagnetic  Elements 


Element 

Poly/Single  Crystal 

Tc(K) 

Curie  Point  Operation  nRPi 

Ho 

Poly 

20 

50% 

Single 

57% 

Er 

Single 

20 

29% 

Dy 

Poly 

89 

4% 

Single 

23.32% 

Tb 

Poly 

220 

8.5% 

Single 

12.5% 

Gd 

Poly 

293 

10.34% 

Single 

20.5 

Ni 

Poly 

531 

2.5% 

Fe 

Poly 

1044 

9.5% 

Co 

Poly 

1390 

1% 

Y2Fe17 

Poly 

Variable 

9.87% 

An  operating  temperature  differential  (AT)  of  5K  was  considered  about  each  materials  Tc  and  was  used 
to  compute  the  relative  efficiencies;  the  specifics  of  which  can  be  found  in  Hsu  et  al8.  From  the  table  it 
can  be  observed  that  Holmium  (Ho)  and  Erbium  (Er)  maintain  the  largest  relative  efficiencies  when 
oscillating  about  their  Curie  Points.  While  different  reasons  have  been  postulated,  UCLA  believes  it  is 
the  transition  from  a  ferromagnetic  to  antiferromagnetic  state,  using  an  ordered  to  ordered  transition, 
rather  than  a  typical  ferromagnetic  to  paramagnetic  transition  (i.e.  order  to  disorder).  The  transition  of 


materials  from  ferromagnetic  to  paramagnetic  results  in  significantly  more  entropic  losses  in  the  system, 
in  turn  causing  a  reduction  of  the  energy  that  can  be  potentially  harvested.  UCLA  began  to  explore  other 
transition  mechanisms  (i.e.  order  to  order)  beyond  those  previously  studied  by  the  classical 
thermomagnetic  community. 

One  such  order-to-order  transition  is  the  Spin  Reorientation  (SR),  a  transformation  which  involves  a 
thermally  induced  rotation  of  the  magnetic  easy  axis  from  one  crystalline  direction  to  another.  Horner 
and  Varma14  first  visited  the  idea  in  1968  and  theorized  that  certain  materials  could  make  a  spontaneous 
transition  from  an  in-plane  to  an  out  of  plane  magnetization  state  due  to  the  temperature  dependence 
of  magnetocrystalline  anisotropy.  Magnetocrystalline  anisotropy  is  defined  as  the  preferential 
alignment  of  a  material's  magnetic  moment  along  a  specific  crystallographic  axis,  which  arises  from 
changes  in  the  crystal-electric  field  as  the  lattice  expands  thermally. 
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Figure  8:  Magnetic  easy  axis  of  gadolinium  as  a  function  of  temperature 


Figure  8,  a  plot  by  Graham15  augmented  to  provide  a  visualization  of  hexagonal  crystal  axes,  shows 
Gadolinium's  magnetization  orientation  (defined  by  theta  (9),  the  deviation  from  the  C-axis)  as  a 
function  of  temperature  for  a  single  crystal.  At  deep  cryogenic  temperatures  the  magnetic  easy  axis  is 
approximately  at  33  degrees  with  respect  to  the  C-axis,  or  approximately  half  way  between  the  basal 
plane  (B-plane)  and  the  C-axis.  As  the  temperature  is  increased  to  approximately  160  C,  the  B-plane 
becomes  the  preferred  orientation.  As  the  temperature  changes  from  230  K  to  250K,  a  noticeable 
magnetic  reorientation  is  observed;  that  is,  the  easy  magnetization  rotates  90°  from  the  basal  plane  to 
the  C-axis.  This  experimental  data  indicates  that  Gd  has  a  Spin  Reorientation  Temperature  (TSr)  at  250° 
K.  If  one  were  to  imagine  that  a  single  domain  element  of  gadolinium  were  wrapped  in  a  theoretical  coil, 
then  as  the  element  was  heated  from  230  to  250K,  the  magnetetization  would  rotate  from  the  B  axis  to 


the  C  axis  and  this  flux  change  would  generate  current  in  the  wire,  i.e.  electrical  energy.  Utilizing  this 
foundation,  we  began  exploring  the  integration  the  SR  phenomenon  into  a  harvesting  methodology. 


3.0  Spin  Reorientation:  Bulk  Modeling 


Our  first  task  was  to  model  the  thermal  to  magnetic  transduction  mechanism,  which  was  initiated  by 
modeling  a  system  at  the  bulk  scale.  We  chose  single  crystal  Gd  as  the  working  element  and  integrated 
its  magnetocryalline  parameters,  found  experimentally  by  Graham15 ,  into  the  COMSOL  finite  element 
package  magnetics  module.  The  ideal  geometry  for  a  flux  yoke,  given  the  orthogonal  magnetization 
paths  during  a  spin  reorientation  transition,  was  provided  by  Ohkoshi16 . 
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Figure  9:  Left  -  The  proposed  magnetic  flux  return  yoke  by  Ohkoshi;  Right  -  The  system  modeled  in  COMSOL  above  TSR 


Figure  9  Left,  illustrates  Ohkoshi's  original  design  showing  that  below  the  SR  temperature  the  magnetic 
flux  would  flow  along  the  right  loop  of  the  "pretzel",  corresponding  to  the  B-plane  direction  of  the 
NdCo5  crystal,  and  would  switch  to  the  left  path,  corresponding  to  the  C-axis,  as  the  temperature  was 
increased  beyond  the  SR  point.  Coils  1  and  2,  wrapped  around  the  left  and  right  paths  respectively, 
would  see  an  induced  current  as  the  spin  reorientation  occurred.  In  our  model  we  replaced  the  NdCo5 
with  Gd,  due  to  our  familiarity  and  the  availability  of  material  properties,  and  the  SmCo5  biasing  magnet 
with  NdFeB,  a  hard  magnet  which  provides  a  larger  magnetic  field  which  was  unavailable  during 
Ohkoshi's  study.  In  a  system  such  as  this,  the  amount  of  voltage  induced  within  a  theoretical  coil, 
wrapped  around  each  path,  is  monotonically  related  to  the  permeability  change  of  the  material  which 
makes  up  the  flux  path  by  Lenz's  Law. 


In  this  equation  p  is  the  permeability,  t  is  time,  A  is  the  cross-sectional  area  of  the  flux  path,  H  is  the 
magnetic  field  and  L  is  the  length  of  the  magnetic  circuit.  As  the  area,  path  length  and  applied  field 
remain  invariant,  only  the  permeability's  change  contributes  to  the  energy  output.  Magnetic  steel  was 
chosen  as  flux  path  material  in  this  study  due  to  its  relatively  high  permeability  under  the  applied  bias 
field  of  3  kOe,  such  that  the  spin  reorientation  in  the  Gd  would  induce  the  largest  possible  change  in  the 
yoke. 


Figure  10:  Permeability  of  the  steel  yoke  as  a  function  of  applied  field  at  217K  along  both  flux  paths 


Figure  11:  Permeability  of  the  steel  yoke  as  a  function  of  applied  field  at  237K  along  both  flux  paths 

The  permeability  of  the  yoke  in  each  path  direction  was  investigated  for  each  temperature  in  the  SR 
region.  Figures  10  and  11  show  the  permeability  of  the  yoke  in  the  spin  reorientation  zone  at  217  and 
237K  respectively.  It  should  be  recognized  that  in  the  field  region  required  to  saturate  Gd  in  either 
direction,  highlighted  by  the  blue  box,  that  the  permeability  of  each  flux  path  changes  the  most  across 
the  SR  temperature  differential,  flipping  from  B-plane  dominated  to  C-axis.  This  region  was  chosen  as 
the  operating  zone  for  further  study. 


Figure  12:  The  permeability  change  of  the  yoke  across  gadoliniums  SR  region  at  an  applied  field  of  H=3200  Oe 


Figure  12  shows  the  change  in  permeability  as  a  function  of  temperature  at  the  saturating  field.  From 
the  figure  it  can  be  seen  that  as  the  magnetization  of  the  gadolinium  rotations  90°,  that  the  relative 
permeability  of  the  steel  yoke  changes  by  approximately  10.  Using  Lenz's  law  to  back  out  an  induced 
voltage  yields  an  energy  density  of  5.1  kJ/m3.  With  this  value  being  only  a  fifth  of  the  energy  predicted 
when  operating  about  Tc,  we  took  a  step  back  to  analyze  the  thermodynamics  of  the  spin  reorientation 
transition. 

4.0  Thermodynamic  Modeling  of  Spin  Reorientation 

In  order  to  determine  the  ideal  parameters  for  harvesting  energy  using  a  SR  transition  we  required  a 
better  understanding  of  the  thermodynamics  that  occur  as  magnets  are  thermally  cycled.  A  model  was 
developed  to  describe  how  each  magnetic  contribution  to  the  total  energy  varied  with  temperature  with 
considerations  given  to  the  different  thermodynamic  processes  that  occur  across  Tc  and  SR.  The 
following  section  outlines  the  formulation  used  to  create  this  model. 

In  the  most  basic  magnetic  system,  the  magnetic  contribution  of  the  energy  is  E0  =  )  where  HT 

is  the  total  field  and  p0  is  the  permeability.  If  one  assumes  that  at  an  initial  state  there  is  no  magnetic 
medium,  then  the  total  magnetic  energy  in  a  system  results  from  an  externally  applied  field  (HA)  only.  If 
the  system  is  modified  to  include  a  magnetic  material  and  the  applied  field  is  considered  to  be  external 
to  the  system,  then  the  magnetic  energy  of  the  enclosed  system  can  be  represented  as 

E  =  -H^y  Here  the  difference  between  the  squares  of  the  field  quantities  represents  the 

local  energy  encapsulated  in  a  control  surface  surrounding  the  boundary  of  the  magnetic  media.  That 
field  quantity  HT  has  contributions  due  to  the  applied  external  field  and  a  demagnetization  field  (HD) 
which  acts  to  oppose  the  external  field  within  the  magnetic  material,  resulting  in  the  following. 

E  -  ^((#,  -  ■ Hd  f  -  H\ ) -  (X ' + Hi  +  2 HaHd  -  X)  (4.D 


After  expanding  and  simplifying  the  result,  the  relation 


hd- 


is  applied  to  obtain  the 


\f*o  J 

"demagnetization  energy"  in  terms  of  the  magnetization  (M)  of  the  magnetic  material  and  flux  due  to 
the  demagnetization  field  (BD).  This  formula  is  then  expended  as  follows. 
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The  first  term  of  the  expanded  Equation  4.2  can  be  represented  as  the  following  volume  integral  in  all 
space: 


,  ft 

Where  the  demagnetization  field  HD  is  equal 
Equation  4.5. 
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to  the  negative  gradient  of  a  scalar  field,  -V0,  to  yield 
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(4.4) 


By  using  the  electrostatic  Maxwell  relation,  that  the  divergence  of  B  is  equal  to  zero,  the  second  term 
can  be  reduced  to  zero.  From  there,  the  volume  integral  over  all  space  can  be  converted  into  a  surface 
integral,  by  Gauss'  Law,  to  yield  Equation  4.6. 


—  f((^BD)g,)dS 

Vo  A.S. 

As  the  surface  approaches  infinity,  such  that  lim(r3^)  — >  0  and 

r—>cc  '  t 

(f)D  -*  0  ,  the  entirety  of  Equation  4.3  can  be  seen  to  approach  0. 
to  yield  Equation  4.6. 
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j  — *  0,  which  implies  that 
4.2  can  then  be  simplified 
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The  relation  H D  = 
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is  then  applied  to  the  second  term,  in  order  to  have  both  energies  in 


--M 

\V0  , 

terms  of  magnetization.  Simplifying  the  total  energy  transforms  Equation  4.6  into  the  following. 


E  =  - 


+  HaBd-/u0HaM 
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(4.7) 


Once  again  the  product  of  a  field  and  a  flux,  in  this  case  HABD,  goes  to  zero  by  invoking  electrostatic 
Maxwell  relations  and  Gauss  theorem.  Then  equating  the  demagnetization  field  (HD)  to  the  product  of  a 

shape  factor  matrix  (N)  and  the  magnetization  such  that HD  =  —NM ,  provides  a  continuum  energy 
function  for  a  magnetic  system  as  Equation  4.8. 

E  =  ^-M(NM)-voHaM  (4.8) 

From  the  Hamiltonian  describing  the  magnetic  system  we  know  that  two  additional  energy  terms  should 
be  added  to  describe  the  quantum  effects  acting  on  magnetic  system;  one  to  describe  the  exchange 
energy  (Eex)  and  one  to  describe  the  magnetocrystalline  anisotropy  energy  (EMca)  arising  from  spin-orbit 
interactions  and  the  crystal-electric-field. 

E.^-M(NM)-^HaM+Eei+Exa  (4.9) 

Under  the  assumption  of  a  single  domain  magnetic  element,  we  can  consider  any  rotation  of  the  spins 
to  be  coherent  and  thus  any  two  states  will  have  identical  exchange  energies  and  we  will  henceforth 
drop  the  term.  The  magnetocrystalline  anisotropy  energy  can  take  multiple  forms  depending  of  the 
point  group  of  the  crystal  being  considered.  In  this  study  we  evaluated  Gd  and  NdCo5,  both  HCP  crystals 
who  have  EMca  of  the  form: 

Emca  =AT1sin2(^)  +  ir2sin4((9)  +  ^r3sin6((9)Ar4sin6(^)cos(6^)  (4.10) 

Using  our  single  domain  assumption  and  applying  Equation  4.10  for  the  magnetocryalline  anisotropy 
yields  a  total  magnetic  energy  as  follows. 

Em  =  ~M  ( NM )-  /u0HaM  +  sin2  (d)  +  K2  sin4  (d)  +  K2  sin6  ( d)K4  sin6  (tf  )cos(6^)  (4.11) 

From  this  equation  it  can  be  seen  that  the  terms  can  be  described  with  independent  magnetization  [M] 
and  temperature  [T]  which  indicates  that  the  magnetic  energy  is  a  free  energy  of  the  Helmholtz  variety: 

F  =  U  -TS  =  Em  - SdT .  One  must  then  add  to  that  the  energy  associated  with  the  crystal  lattice.  Due 

to  considering  temperature  as  an  independent  variable,  one  must  account  for  the  effects  of  thermal 
energy  on  the  lattice,  via  the  temperature  dependent  heat  capacity  of  the  system. 


dQ  =  CpdT^Q=foCp{T)dT 


(4.12) 


In  the  transition  from  an  internal  energy  to  a  free  energy  this  term  takes  on  a  negative  sign  from  the 
standard  Legendre  transform  yielding  the  Helmholtz  free  energy  as  Equation  4.12  subtracted  from 
Equation  4.11.  For  ease  of  measurement  we  want  a  thermodynamic  system  in  independent 
temperature  [T]  and  magnetic  field  [H],  the  variables  most  easy  to  prescribe,  we  transform  the  system 

into  a  Gibbs  free  energy  by  the  following:  G  =  F  +  / lQHAM .  This  yields  a  total  Gibbs  free  energy  as 

G  =  —M(NM)  +  K,  sin2  (ff)  +  K2  sin4  (o)  +  K3  sin6  sin6  (tf  )cos(6^)  -f  C  ( T)dT  (4.13) 
2  4° 

In  order  to  quantify  the  work  done  per  cycle,  the  total  entropy  of  the  system  is  defined  as  the  negative 
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partial  derivative  of  the  Gibbs  free  energy  with  respect  to  temperature,  S  = - - - - . 
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(4.14) 


The  total  energy  change  is  thus  equivalent  to  the  amount  of  heat  that  flows  into  the  system:  Q  =JSdT 
,  where  the  magnetic  change  under  the  variation  in  temperature  constitutes  the  work:  W  =  JSmagdT. 


Note  that  the  magnetic  contribution  of  the  entropy  includes  every  term  except  for  the  lattice  heat 
capacity  Cp.  The  absolute  efficiency  of  the  system  can  then  be  denoted  as  the  magnetic  work  over  the 
total  heat  input. 
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As  thermodynamic  cycles  are  subject  to  the  Carnot  limit,  the  ratio  of  the  high  temperature  over  the 
change  in  temperature,  heat  engines  are  typically  compared  by  their  efficiency  relative  to  that  limit  by 
the  following. 


Vrel  = 


(4.16) 


Using  this  standardized  metric  as  well  as  the  total  magnetic  energy  density  computed  by  Equation  4.15, 
Gd  and  NdCo5,  operating  in  their  spin  reorientation  regimes,  were  evaluated  against  one  another  as  well 
as  Gd  operating  about  its  Curie  point. 

5.0  Analytical  Results 

In  order  to  compute  the  energy  outputs  and  efficiencies  of  the  relevant  systems,  Gd  and  NdCo5,  the 
thermodynamic  energy  formulation  was  coded  into  Mathematica  to  facilitate  a  numerical  solution  for 
each  process  in  the  thermodynamic  cycle.  This  was  accomplished  by  first  obtaining  the  saturation 
magnetization  (Ms),  heat  capacity  (Cp),  and  magnetocrystalline  anisotropy  coefficients  (K's)  as  functions 
of  the  temperature  from  the  literature11'14-16  and  fitting  each  set  with  a  spline,  for  continuous 
differentiability.  Next  we  assumed  that  the  shape  of  our  magnetic  elements  were  ellipsoidal,  and  thus 
applied  Osborne  shape  factors14  into  the  demagnetization  tensor  (N).  With  values  for  every  material 
and  geometrical  parameter  coded  in,  we  then  went  about  computing  the  location  of  the  magnetization 
as  a  function  of  temperature  and  applied  field. 


Figure  13:  Left  -  Crystal  orientation  showing  the  applied  field  vector  (H)  in  green  and  its  corresponding  polar  (a)  and 
azimuthal  (p)  angle  as  well  as  the  red  magnetization  vector  (M)  and  its  polar  (6)  and  azimuthal  (4>)  angles.  In  addition  the 
black  and  blue  arrows  indicate  the  C  and  B  axes  respectively.  Middle  -  The  energy  surface  of  a  NdCo5  sphere  at  T=255K  with 
a  100  Oe  field  applied  at  a  =  45°  and  p  =  45°.  Right  -  The  energy  surface  of  NdCo5  at  286K  with  the  same  applied  field 
parameters. 

The  numerical  energy  minimization  used  to  determine  the  magnetization  direction  was  computed  from 
the  total  magnetic  energy,  Equation  1.12,  by  taking  its  partial  derivative  relative  to  the  magnetization. 


Due  to  the  hexagonal  crystal's  six  fold  symmetry,  this  global  minimization  had  the  potential  to  return  six 
magnetization  vectors  which  were  all  equivalent  in  the  absence  of  symmetry  breaking  applied  field  or 
shape  effects.  It  was  possible  to  obtain  a  set  of  quasi-static  magnetization  results  which  were 
discontinuous  across  a  set  of  temperature  values.  Figure  13  provides  two  examples  of  the  energy 
surfaces  for  NdCo5.  It  is  apparent  that  NdCo5  has  an  energy  surface  with  6  symmetrical  energy  wells, 
noting  that  blue  surface  corresponds  to  a  negative  total  energy  value  and  red  to  positive  values; 
negative  energy  values  are  a  result  of  the  magnetocryalline  parameters,  which  are  not  normalized  to  0 
at  absolute  zero  but  rather  at  Tc  when  magnetic  cohesion  disappears.  According  to  experimental 
studies15,  the  magnetization  vector  (shown  in  red)  should  move  continuously  from  one  set  of  angles  to 
another.  To  remedy  the  potential  discontinuities  in  the  numerics,  we  broke  the  energy  minimization  into 
two  parts. 


Figure  14:  Magnetization  angle  for  NdCo5  through  its  spin-reorientation  zone 


The  first  part  was  comprised  of  the  global  energy  minimization,  using  Mathematica's  built  in 
"RandomSearch"  method,  to  obtain  a  global  minima  at  the  starting  temperature  of  interest.  The  second 
part  started  from  the  solution  of  the  first  and  traversed  the  energy  landscape  to  the  nearest  equivalent 
minima  at  the  subsequent  temperature  values.  This  methodology  guaranteed  a  continuous  table  of 
magnetization  vectors,  the  axial  angle  is  plotted  in  Figure  14  for  NdCo5.  Below  the  SR  region,  the  6 
energy  minima  lay  in  the  basal  plane  and  as  the  temperature  increases  beyond  260K  two  easy  cones 


form,  while  energy  maxima  form  along  B-plane.  As  the  temperature  is  increased  beyond  the  SR  region 
the  energy  minima  collapses  to  the  C-axis.  In  simple  language,  the  easy  axis  rotates  from  90  degree 
below  260  K  to  0  degree  above  286  K.  With  continuous  functions  of  the  magnetization  direction 
computed,  it  was  then  possible  to  calculate  the  thermomagnetic  work  potential  and  efficiency  of  this 
transition. 


Table  2:  Energy  and  Efficiency  Results 


Material 

Transition 

Tl  [K] 

Th  [K] 

Qin  [J/m3] 

Wout  [J/m3] 

Habs 

Hrel 

Gd 

SR 

219 

243 

4.67  E7 

10.2  E5 

0.02% 

0.01% 

Gd 

Curie 

283 

288 

1.31  E7 

6.83  E5 

0.52% 

30.4% 

NdCo5 

SR 

255 

286 

1.04  E8 

2.57  E6 

2.4% 

22.2% 

Using  Equations  1.13  and  1.15  through  1.17  and  the  numerically  computed  magnetization  states,  the 
work  and  efficiencies  were  compiled  into  Table  2  for  Gd  spin  reorientation,  Gd  Curie  point  operation, 
and  NdCo  spin  reorientation.  In  these  calculations  the  applied  magnetic  field  was  100  Oe  at  a  and  (3 
equal  to  45°,  this  guaranteed  that  the  spin  reorientations  always  occurred  in  similar  directions. 
Additionally  the  operating  temperature  differential  was  chosen  to  increase  the  relative  efficiency  (nrei) 
for  the  given  calculation  rather  than  keep  the  temperature  differential  constant  for  each  study.  For  Gd  it 
can  be  seen  that  operating  about  its  SR  temperature  actually  provides  a  significant  reduction  in 
performance  compared  to  its  Curie  point;  compare  30%  relative  efficiency  for  Curie  to  0.01%  relative 
efficiency  for  spin  reorientation.  NdCo5,  on  the  other  hand,  showed  an  order  of  magnitude  increase  in 
work  density  increase,  as  well  as  an  increase  in  absolute  efficiency  as  compared  to  either  Gd  example  in 
the  Table  2.  The  relative  efficiency,  however,  was  about  8%  smaller  than  the  ideal  Gd  Curie  point 
transition;  this  appears  reasonable  as  the  operating  temperature  differential  is  approximately  6  times 
smaller  for  Gd  (i.e.  5  K)  about  its  Tc  than  NdCo5  (i.e.  31  K)  about  its  SR  which  significantly  reduces  the 
Carnot  limit.  We  checked  the  validity  of  the  calculations  by  determining  which  term  dominated  Equation 
1.15  across  the  operating  temperature  differential  and  compared  its  magnitude  change  to  the 
experimental  results  available  in  the  literature. 


Figure  15:  Left  -  Magnetocryastalline  paramerters  of  Gd  from  Graham10  .  Right  -  Magnetocrystalline  parameters  of  NdCo5 
from  Klein16 .  Each  is  demarked  with  the  approximate  energy  change  in  the  SR  zone. 

When  the  materials  operate  about  their  SR  zones,  the  magnitude  of  their  magnetization  is  nearly 
constant,  and  the  majority  of  the  temperature  induced  change  takes  place  in  the  magnetocrystalline 
constants.  Figure  9  provides  the  values  for  these  parameters  as  functions  of  temperature.  The 
temperature  dependent  magnetocrystalline  change  is  approximately  on  the  same  order  as  the 
calculated  work.  For  operation  about  the  Curie  point  the  magnetocrystalline  terms  reduce  to  nearly 
zero,  such  that  only  the  magnitude  of  magnetization  changes  with  temperature.  The  results  for  Gd 
about  its  Curie  point,  computing  work  from  entropy,  match  the  results  in  Hsu  et  al8  where  the  work  was 
computed  from  M-H  loops. 

A  note  should  be  made  about  the  difference  between  the  calculated  results  for  NdCo5  in  comparison  to 
results  previously  reporting  a  relative  efficiency  approaching  60%.  In  our  original  computation  an  error 
was  made  in  determining  the  required  heat  input  from  the  specific  heat  (Cp).  This  term  was  reported  by 
Ohkoshi16  in  units  erg/(mol  K)  which  was  inappropriately  converted  to  erg/(g  K),  instead  of  erg/(cm3  K), 
due  to  lack  of  information  regarding  the  density  of  NdCo5.  A  density  of  8.35  g/cm3  was  determined  by 
computing  a  hexagonal  unit  cell  from  the  bond  lengths  and  atomic  sizes  that  were  provided.  This  value  is 
within  10%  of  computing  the  density  by  a  rule  of  mixtures  between  the  two  elements  and  is  thus 
considered  relatively  accurate.  We  postulated  that  even  with  the  reduction  in  relative  efficiency,  that 
the  order  of  magnitude  increase  in  work  and  absolute  efficiency  warranted  attempting  to  fabricate 
nano-scale  elements  to  experimentally  verify  our  calculations.  The  amount  of  energy  being  produced 


and  the  "decent"  relative  efficiency  of  NdCo5  still  makes  it  an  attractive  candidate  for  thermal  energy 
harvesting  based  on  this  work. 

6.0  Spin  Reorientation  Experimental  Effort 

The  following  paragraphs  focus  on  sputter  deposition  of  thin  film  epitaxial  NdCo5.  While  this  was 
attempted  we  did  not  achieve  our  stated  objective  of  creating  an  epitaxial  thin  film  of  NdCo5.  Unlike 
elemental  gadolinium,  NdCo5  is  a  bimetallic  compound  which  requires  precise  stoichiometric  control  to 
achieve  the  desired  magnetic  properties.  In  the  literature  epitaxial  thin  films  have  been  produced  by 
pulsed  laser  deposition  (PLD)  using  co-deposited  single  element  targets18,  however  this  method  is 
currently  unavailable  to  UCLA  and  in  the  potential  event  of  desired  scale-up  would  prove  to  be  cost  and 
time  prohibitive.  Thus  we  attempted  to  epitaxially  deposit  this  material  by  co-sputtering. 
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Figure  16:  EDS  results  of  Nd-Co  deposition  on  Si  substrate.  Inset  shows  co-sputtering  configuration. 

Neodymium  and  cobalt  targets  of  99.9%  purity  were  obtained  from  the  Kurt  J.  Lesker  Company  and 
calibrated  into  our  Denton  UHV  sputtering  system.  Figure  12  shows  the  co-sputtering  configuration  and 
our  initial  attempt  to  deposit  the  material  on  a  silicon  substrate.  The  parameters  of  deposition  were  an 
Ar  pressure  and  flow  rate  of  2.5  mTorr  at  20  SCCM,  sputtering  powers  of  280  and  42  Watts  for  Nd  and 
Co  respectively,  a  pre-sputter  time  of  360  sec,  and  a  target  to  substrate  distance  of  5cm.  After 
depositing  the  100  nm  film,  the  composition  was  examined  using  electron  diffraction  spectroscopy  (EDS) 


to  provide  the  atomic  percentage  of  each  element  present.  Stoichiometric  NdCo5  should  be  16.67%  Nd 
and  83.33%  Co,  but  our  initial  attempt  showed  approximately  5%  more  Nd  content  than  desired. 
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Figure  17:  XRD  results  of  recrystallization  study  on  original  Nd-Co  deposition 

Even  with  this  disparity  we  performed  measurements  to  create  a  baseline  for  recrystallization 
temperatures.  This  work  included  measuring  the  XRD  response  following  annealing  the  sample  in 
atmosphere  at  temperatures  ranging  between  500  and  850  C.  The  results  are  presented  in  Figure  13. 
The  data  shows  that  the  sample  became  crystalline  once  annealed  above  750  C.  While  a  moderate  peek 
presented  itself  along  the  (002)  C-axis  direction,  the  largest  peak  was  in  a  cubic  (101)  direction.  This  can 
be  partially  attributed  to  compositional  variation,  but  we  suspected  that  the  cubic  Si  substrate  provided 
the  greatest  contributing  factor  and  one  that  would  not  produce  epitaxial  NdCo5  structures.  However, 
these  results  indicate  that  temperatures  above  750  C  are  required  to  crystallize  the  deposited  films. 


Figure  18:  Left  -  EDS  of  100  nm  of  NdCo5  film  on  sapphire.  Right  -  XRD  of  100  nm  NdCo5  film  on  sapphire. 

In  subsequent  runs  we  used  a  sapphire  substrate,  a  hexagonal  crystal  with  lattice  parameters  within  4% 
of  the  desired  NdCo5  crystallography,  to  deposit  the  NdCo5  film.  This  substrate  should  produce  an 
epitaxial  film  if  the  deposition  rate  and  crystallization  rate  is  adequately  be  controlled.  During 
deposition  on  this  substrate  we  conducted  a  parametric  study  proportionally  modifying  the  driving 
power  to  the  plasma  sources  to  tune  the  NdCo  composition.  Figure  14  provides  the  EDS  and  XRD  results 
for  a  NdCo5  sample  once  the  appropriate  deposition  parameters  were.  From  the  EDS  results,  large  peaks 
appeared  corresponding  to  the  aluminum  and  oxygen  present  in  the  substrate,  however  it  can  be  seen 
that  after  subtracting  out  the  counts  due  to  the  substrate  that  the  desired  NdCo5  composition  was 
obtained  (i.e.  see  circled  composition  in  Figure  14).  The  XRD  study  in  Figure  14  right  was  conducted  on  a 
both  bare  sapphire  wafer  and  the  unannealed  film  sample  which  showed  very  strong  agreement 
between  the  lattice  structures  of  the  two  materials,  indicating  that  film  was  apparently  epitaxial 
deposited  without  post  deposition  annealing.  This  result  suggested  that  the  deposition  energy  was 
sufficient  to  at  least  partially  crystalize  the  film  as  a  crystalline  structure  is  observed  in  the  XRD. 


Figure  19:  Magnetization  curves  of  100  nm  of  NdCo5  on  sapphire 

Using  the  unannealed  sample  that  appeared  to  be  epitaxial  deposited  shown  in  Figure  14,  the  NdCo5 
film  magnetic  response  was  measured  using  SQUID  at  temperatures  below  (240)  and  above  (298  and 
350)  the  SR  zone.  The  results  are  presented  in  Figure  15.  The  magnetization  curves  were  taken  in  the 
out-of-plane  (C-axis)  direction.  It  was  expected  that  the  M-H  loops  would  show  a  hard  axis  at  less  than 
250  K,  which  would  then  transition  to  an  easier  axis,  with  a  larger  coercivity,  above  280  K.  Here  it  is 
important  to  also  point  out  that  the  overall  coercivity  was  only  200  Oe,  much  lower  than  the  expected 
2500  Oe  at  350  K  in  a  single  crystal.  We  originally  attributed  the  smaller  coercive  field  to  the  lack  of  a 
fully  crystalline  material;  we  hypothesized  that  the  materials  contained  both  amorphous  and  crystalline 
regions  through  the  thickness  of  the  film  producing  a  softer  magnetic  structure.  To  resolve  the 
amorphous/crystalline  problem  we  annealed  a  portion  of  the  sample  at  850  C.  However,  upon  magnetic 
measuring  of  the  annealed  sample  and  subsequent  measurements  on  composition,  we  found  that  the 
Nd  content  of  the  annealed  portion  dropped  substantially.  The  remaining  thin  film  contained  only  trace 
Nd  content.  At  this  stage  we  determined  that,  like  Gd,  a  capping  layer  would  be  necessary  to  prevent 
loss  of  Nd. 


Figure  20:  Magnetization  curves  for  100  nm  of  NdCo5  on  sapphire  with  10  nm  Ti  capping  layer 

During  the  next  set  of  depositions  10  nm  of  Ti  was  sputtered  over  the  top  of  the  NdCo5  in-situ,  before 
breaking  vacuum.  Ti  was  chosen  because  this  target  was  in  the  shared  resource  facility  already  present 
in  the  sputtering  system.  The  argon  flow  rate  was  reduced  from  20  SCCM  to  12  SCCM  in  an  attempt  to 
slow  the  deposition  rate  and  attain  a  higher  ordering  of  the  NdCo5  crystals.  Figure  16  shows  the 
magnetization  curves  for  the  Ti  capped  sample.  Using  the  same  experimental  procedure  in  the  SQUID 
with  additional  measurements  performed  at  colder  temperatures,  a  spin  reorientation  appears  to  be 
present  between  200  and  350  K.  While  this  temperature  differential  is  much  larger  than  expected,  it 
showed  that  this  deposition  method  had  promise  to  achieve  the  desired  effect.  Additionally,  we 
hypothesized  that  the  lower  coercivity  value,  100  Oe  compared  to  the  predicted  2500  Oe,  was  due  to 
internal  stresses  and  could  be  rectified  by  annealing  the  sample  to  relax  the  crystal  lattice. 


Figure  21:  Left  -  Magnetization  curves  for  100  nm  NdCo5  on  sapphire  with  10  nm  Ti  and  annealed  to  850  C 


This  sample  was  subsequently  annealed  to  850  C,  for  1  hour,  in  a  tube  furnace  under  a  reducing 
environment,  90%  argon/10%  hydrogen  to  prevent  oxidation.  Figure  17  shows  the  results  of  the 
repeated  magnetization  curves  after  annealing  and  subsequent  XRD.  It  can  be  seen  that  the  annealing 
increased  the  coercivity  from  approximately  80  to  nearly  250  Oe,  which  was  attributed  to  a  reduction  in 
film  stress  and  potential  relaxation  of  the  crystals  into  a  more  ordered  state.  Unfortunately,  the 
reorientation  observed  before  annealing  was  absent.  The  XRD  appears  to  indicate  that  post  annealing 
relaxes  a  portion  of  deposited  film  into  a  cubic  form,  with  a  notable  peak  forming  in  the  (110)  direction 
as  well  as  the  hexagonal  (002).  Additionally,  WDS  results  showed  that  the  Ti  capping  layer  was  depleting 
the  Nd  during  the  annealing  process.  To  prevent  this  reaction  in  subsequent  depositions,  future  runs 
were  capped  with  tantalum  (Ta),  a  material  which  is  much  more  stable  and  less  reactive  at  high 
temperatures. 
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Figure  22:  Left  -  XRD  results  for  Ta  capped  NdCo5  after  2  hour  anneal  at  1000°  C  for  2  hours;  Right  -  ZFC-FC  results  for  the 
same  sample  between  100  and  350  K. 


Figure  18  shows  the  crystallography  and  temperature  dependent  magnetization  results  for  one  of  the 
stoichiometric  NdCo5  samples  capped  in  Ta  after  a  2  hour  anneal  at  1000°  C.  The  XRD  plot  shows  peak 
counts  very  similar  to  the  Ti  coated  sample  in  Figure  17,  indicating  that  the  sample  crystallized  during 
the  annealing  process.  A  WDS  stoichiometry  study  (not  shown)  indicated  NdCo5  composition  had  been 
achieved.  That  is  the  values  are  close  to  those  shown  in  Figure  14  left,  both  before  and  after  annealing. 
Given  the  compositional  and  crystallographic  success  on  this  sample,  we  examined  the  samples 
magnetization  as  a  function  of  temperature  using  ZFC-FC  and  found  a  valley  at  250  K  during  the  field 
cooling  portion  of  the  experiment.  This  anomaly  in  the  in-plane  magnetization  measurements  suggests 
that  the  magnetization  is  rotating  into  the  direction  of  measurement  at  this  temperature,  matching  the 
theoretical  prediction.  However,  the  magnitude  of  this  change  was  only  approximately  5  emu/cc,  at  an 


applied  field  of  50  Oe.  We  postulated  that  the  magnitude  was  smaller  than  anticipated  due  to  the  fact 
that  this  film  appeared  to  be  polycrystalline  rather  than  textured.  The  less  than  fully  crystallized  regions 
contributed  less  to  magnetization  rotation  due  to  being  less  thermomagnetically  coupled.  Magnetization 
loops  at  constant  temperatures  were  performed  to  confirm  the  hypothesized  lack  of  global  easy  axis 
rotation. 


Figure  23:  Magnetization  curves  for  Ta  capped  NdCo5  annealed  at  1000°  C  for  2  hours. 

The  results  in  Figure  19  confirmed  our  suspicions  that  while  the  sample  was  noticeably  crystalized,  that 
the  uniformity  of  the  crystallization  was  not  epitaxial.  This  can  be  observed  in  the  small  coercivity  values, 
approximately  100  Oe  instead  of  the  expected  2500  Oe,  and  absence  of  a  hard-axis  curve  for  the 
measurements  at  298  and  350  K.  It  was  hypothesized  that  in  order  to  achieve  the  required  crystal 
epitaxy  through  the  thickness  of  the  film,  to  obtain  spin  reorientation,  it  would  be  necessary  to  heat  the 
substrate  at  or  above  850  during  the  co-sputtering  processes;  an  ability  that  UCLA  does  not  currently 
have.  We  attempted  to  reach  out  to  colleagues  at  Stanford  and  Northeastern  and  are  still  in  the  process 
of  doing  this.  More  recently  we  have  attempted  field  alignment,  an  application  of  3000  Oe  along  the  C- 
axis,  during  two  different  annealing  processes:  vacuum  annealing  at  750  C  for  2  hours,  the  max  that  our 
cryo-pump  could  handle,  and  tube  furnace  annealing  at  1000  C  for  1  hour  under  the  reducing 
environment.  Neither  of  these  methods  significantly  improved  the  crystal  uniformity  (i.e.  epitaxy)  or 
magnetic  properties.  The  bottom  line  is  that  we  believe  that  sputter  deposition  can  be  used  to 
successfully  deposit  epitaxial,  thin  film  NdCo5.  However,  in-situ  substrate  heating  at  850°  C  is  required 
with  both  minimal  Ar  pressure  and  under  UHV.  The  substrate  heating  provides  sufficient  mobility  of  the 
co-sputtered  material  to  nucleate  textured  crystal  growth,  deposit  uniformly  without  large  stresses,  and 
remain  contaminant  and  oxide  free.  While  this  was  not  successfully  achieved  during  this  program,  we 


are  continuing  this  study  in  the  next  couple  of  months  with  vacuum  annealing  to  be  performed  at 
Brookhaven  National  lab. 


Conclusion 

Thermomagnetics,  and  specifically  spin  reorientation,  show  promising  advantages  over  other  thermal 
harvesting  methods,  namely  an  order  of  magnitude  increase  in  work  output  and  absolute  efficiency. 
The  thermodynamic  model  developed  under  this  program  represents  an  advancement  for  the 
community  in  studying  this  fairly  unexplored  area.  The  model  predicts  the  theoretical  work  output  from 
entropy  and  provides  an  analytical  approach  to  compare  different  2nd  order  thermomagnetic 
transformations.  Future  improvements  to  the  model  could  include  discontinuous  1st  order  systems. 
While  the  results  computed  for  NdCo5  are  promising,  when  compared  to  Gd,  they  still  fall  short  of  the 
50%  relative  efficiency  mark  we  had  originally  hoped  to  obtain.  They  do,  however,  confirm  one  of  our 
hypotheses  that  thermomagnetic  transduction  mechanisms  of  the  order-to-order  type  have  the 
potential  to  achieve  a  higher  absolute  efficiency.  This  is  important  because  this  work  shows  that  the 
field  of  thermomagnetics  is  promising  but  the  appropriate  material  has  not  yet  been  found  that 
produces  a  relative  efficiency  sufficiently  larger  than  currently  exists.  Finally,  UHV  Sputtering  shows 
potential  as  a  means  of  depositing  crystalline  thermomagnetic  nano-structures,  however  additional 
work  needs  to  be  conducted  to  perfect  the  consistency  with  which  these  structures  are  produced. 
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